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Bethe-Salpeter equation in Minkowski space with cross-ladder kernel 

V.A. Karmanov a * and J. Carbonell b 

a Lebedev Physical Institute, 

Leninsky Prospekt 53, 119991 Moscow, Russia 

b Laboratoire de Physique Subatomique et Cosmologie, 
53 avenue des Martyrs, 38026 Grenoble, France 

A new method for solving the Bethe-Salpeter equation is developed. It allows to find the Bethe-Salpeter 
amplitudes both in Minkowski and in Euclidean spaces and, as a by product, the light-front wave function. The 
method is valid for any kernel given by irreducible Feynman graphs. Bethe-Salpeter and Light-Front equations 
for scalar particles with ladder + cross-ladder kernel are solved. 



1. INTRODUCTION 

Bethe-Salpeter (BS) equation £Q is an impor- 
tant tool for studying the relativistic bound state 
problem in a field theory framework (see for re- 
view For a bound state with total momen- 
tum p and in case of equal mass particles, it reads 



*(fc,p) = G& 



(12) 



W) 



jiKi^k 1 ', P )$(k' ; P ) , (i) 



where $ is the BS amplitude, K the interaction 
kernel, m the mass of the constituents, k their rel- 

(12) 

ativc momentum and Gq the free propagators 
of the constituent particles: 

G { Q 12) ( Pl k) = G { 1) (p,k)G i 2) (p,k)= (2) 



We will denote by M = 



the total mass of 



the bound state, and by B = 2m — M its binding 
energy. 

The BS equation has singularities which make 
difficult to find its numerical solution. These sin- 
gularities are due to the free propagators J2J but 
can also result from the interaction kernel itself. 
To overcome this difficulty, Wick [3] formulated 
the BS equation in the Euclidean space, by rotat- 
ing the relative energy fco m the complex plane, 

*Supported in part by the RFBR grants 05-02-17482-a 
and 05-02-26615 



i.e., by introducing the variable = —iko. This 
leads to a well defined integral equation - see ltT2"|> 
below - which can be solved by standard meth- 
ods and provides the mass of the system. In this 
procedure the original BS amplitude is however 
lost, and the "rotated" one can no longer be used 
in calculating other physical observables, like for 
instance form factors. 

A successful attempt to calculate the BS am- 
plitude in Minkowski space for the ladder kernel 
was presented in jjj. It was based on the Nakan- 
ishi integral representation of the amplitude 0]. 
Another approach to solve BS equation with sep- 
arable interaction in Minkowski space was devel- 
oped in and applied to the nucleon-nucleon 
system. In jjj, an equation obtained by projecting 
the BS equation on the light-front (LF) plane was 
studied: the LF kernel was found approximately, 
as an expansion in terms of the BS one and the 
original BS amplitude was not reconstructed from 
its LF projection. 

We present in this paper - based on our works 
|8I9| - a new method for solving the BS equation. 
This method provides the BS amplitude both in 
Minkowski and in Euclidean space, depending on 
the value (real or imaginary) of the ko variable. 
Our approach is based on a projection of the BS 
equation on the LF plane, and on the Nakanishi 
integral representation [2J of the BS amplitude. 
The LF projection plays the role of an integral 
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transform which removes the singularities of the 
BS amplitude. The Nakanishi representation re- 
sults in an equation in a closed form and allows 
to easily restore with its solution, the original BS 
amplitude. The transformed equation is derived 
without any approximation and remains equiva- 
lent to the original BS one. Our method is not 
restricted to the ladder kernel. For more com- 
plicated interactions, e.g. cross box, calculations 
become more lengthy, but the additional difficul- 
ties are due to the evaluation of the Feynman di- 
agram for the kernel itself and not to the solution 
of the equation. 

In order to present the method more distinctly, 
we consider the case of zero total angular momen- 
tum and spinless particles. 

In sect. 0we derive the equation allowing to 
find the BS amplitude both in Minkowski and 
Euclidean spaces. In sect. 01 we present corre- 
sponding LF equation and in sect. 21 ~ the BS 
equation in Euclidean space. In sect. [5] the nu- 
merical solutions are presented. Sect. [5] contains 
some concluding remarks. 

2. BETHE-SALPETER EQUATION ON 
THE LIGHT-FRONT PLANE 

Our method is inspired by an existing relation - 
sec cq. (J2J) below - between the (singular) BS am- 
plitude <&(k,p) and the (non-singular) two-body 
LF wave function ip(k±,x). The latter can be ob- 
tained by projecting the BS amplitude on the LF 
plane. We use the covariant formulation ^Hl with 
the LF plane, defined by ui ■ x = with lu 2 = 0. 
The particular choice ui ~ (wo,cD) = (1,0,0,-1) 
results in the standard LF form t + z = 0. We 
apply to both sides of the integral transform: 

/>oo 

d/3$(k + j3w,p) = / df3G ( 12) (k + (3lu, p ) 



<?(tV) 



d 4 k 

W) 



(3) 



In the transformed equation, the BS amplitude 
is written in terms of the Nakanishi integral rep- 
resentation |2], which for zero angular momentum 
reads: 

= -==L / dz' / d<y' (4) 



[y ■ 



i 2 - \M 2 



k 2 — p ■ k z' 



A similar representation exists for non-zero angu- 
lar momentum. This representation is valid for a 
rather wide class of solutions. This leads to the 
following equation, derived in 8 , for the weight 
function 5(7, z): 



f 

Jo 



gin'-, z)dj' 



7' + 7 + z 2 m 2 + (1 — z 2 )k? 



d-f' dz' V (7, z; 7', z')g{-y\ z') , (5) 
J-i 



where V is a kernel given in terms of the BS in- 
teraction kernel K by 

V^,z;i,z')= (6) 

dp G { 1] (p, k + I3lo)G ( P (p, k + puj) 

K(k + Puj,k',p)d i k' 



LO ■ p 



[A/ 2 + p ■ k'z' — 7' — k 2 



Gq are the free propagators J5J, the denomina- 
tor arises from the representation Q, the inte- 
gration over k' results from the BS equation Q 
and the integration over (3 from the LF projection 
© . The bound state mass M enters through the 
parameter k 2 = m 2 — \M 2 . 

The four-vectors in the r.h.s of eq. 10 appear 
in form of scalar products, expressed through the 
equation variables (7, z) defined by: 



(7 + z 2 m 2 ) 



k = -■ 



■ k 



p ■ k 



1 — Z 2 U! ■ p 

z[j + z 2 m 2 + (1 — z 2 )k 2 
l-z 2 



—z 

2 



and related to the standard LF variables by 7 = 
k\, z = 1 - 2x . 

Equation (JjjJ is equivalent to the initial BS 
equation Q and provides, for a given kernel K, 
the same bound state mass M. Once (7(7, z) is 
known, the BS amplitude can be restored by eq. 
(@J. For real values of kg, we find the BS ampli- 
tude in Minkowski space and for imaginary values 
ko = iki, in the Euclidean one. The correspond- 
ing LF wave function tp(k±,x) is easily obtained 



3 



by 



x(l — x) 



(7) 



$(fc + [3u,p)(uj -p)d/3 = 
x(l- x)g(~f',l-2x)d~f' 



7' + k\ + m 2 - x(l - x)M' 

2.1. Ladder kernel 

Application to the ladder BS kernel 

allows to test our method. In this case, the kernel 
V = of equation (JSJ) is calculated analyti- 

cally. 

It becomes especially simple for fi = 0. This 
particular case constitutes the Wick-Cutkosky 
model [311 lj . We search the solution of (JSJ in 
the form: 

3(7,2) = ^(7) g(z) 

justified by possibility to find a function g(z) 
which does not depend on 7. The integration 
over 7' in both sides of equation JSJ) drops out, 
the 7-dependence of the kernel V is reduced to a 
factor which is the same in both sides of equation 
and cancels. After that, the equation |J5J takes a 
simplified form 

dz'V(z,z')g(z') (9) 



a 

,2 




with a = g / (16nm ) and 

V(Z,z) = : 5 

m 2 - \{l - z' 2 )M 2 

if -1 < z' < z < 1 
if -1 < z < z' < 1 

which exactly coincides with the Wick-Cutkosky 
equation [21311 lj . 

Numerical results for non-zero jj, are given in 
sect. |3J 

2.2. Cross ladder kernel 

The cross-ladder BS kernel is shown in Figure 
n The corresponding amplitude reads: 



K^ CL \k,k',p) 



\p" — m 2 + ie] [(p' 2 — pi + p") 2 — m 2 + ie] 
d 4 p" 



[(Pi ~ P") 2 ~ fJ- 2 + «e] Wx - P") 2 - A' 2 + *e] 
We calculate this expression, substitute the re- 
sult in © , calculate the integrals and find in this 
way the cross-ladder contribution to the kernel 
V(7, z; 7', z') in equation ©. The full kernel is 
the sum of ladder and cross-ladder graphs: 

V = V (L) + y{CL) 

The details of calculation and expression for 
y( CL ) are given in HJ. 



P,-P 



p,-p,+p 



Figure 1. Feynman cross graph. 



3. LIGHT-FRONT EQUATION AND 
KERNELS 

We would like to compare the results obtained 
in the BS approach with the ones found in Light- 
Front Dynamics (LFD). With this aim, we pre- 
cise here the LF equation and the corresponding 
kernel. In terms of k± and x variables, the LF 
equation reads (see e.g. [rU]): 

I k 2 , +m 2 9 \ . -> , m 2 . , 



x / V LF (k' ± ,x';k ± ,x,M 2 ) 



(2tt)< 



(10) 



2jj(k' ±1 x')d 2 k' ± dx' 
2x'(l - x 1 ) ' 

There are two time-ordered ladder graphs and 
six LF time-ordered cross-ladder graphs. The lat- 
ters have the order a 2 . In addition, and to the 
same order a 2 , there are two irreducible time- 
ordered graphs with two mesons in the intermedi- 
ate state (stretched boxes). The full LFD kernel - 



4 



including ladder, cross-ladder and stretched-box 
graphs - is written in the form: 



V LF (k' ± , x'; k ± ,x, M 2 ) = V$ + Vt F L) + V, 



LF 



The LF equation l|llf) has been solved with this 
full kernel. 



5. NUMERICAL RESULTS 

It turns out that the discretized integral opera- 
tor in l.h.-side of eq. 10) has very small eigenval- 
ues. They are unphysical but make the solution 
unstable. We have regularize it by adding a small 
constant e to its diagonal part. This procedure 
allows us to obtain stable eigenvalues in the in- 
terval e = 1CT 4 1CT 12 . 



4. CROSS LADDER IN EUCLIDEAN 
SPACE 

The possibility of Wick rotation has been 
proved for the ladder kernel. One can also " ro- 
tate", without crossing singularities, simultane- 
ously all the energies kio of a perturbative Feyn- 
man amplitude. However, this possibility is not 
evident for a BS amplitude with a kernel more 
complicated than the ladder one, since now not 
all the energies are rotated: when the relative en- 
ergy k is replaced by ik 4 , the eigenvalue - the 
bound state energy pq - still remains real. For 
the latter case, the validity /invalidity of Wick 
rotation is discussed in the literature. For the 
cross-ladder kernel, we can check it numerically, 
by finding the binding energies from eq. (0 - 
see results on Table |5] below - and by solving 
the Euclidean space equation for L+CL kernel. 
One can then see whether or not the Euclidean 
and Minkowski results coincide with each other. 
Coincidence would confirm both: the validity of 
Wick rotation for the cross ladder kernel and our 
method in Minkowski space. 

Equation in Euclidean space is obtained from 
([T]) by the replacement: fco = ik 4 , k' Q — ik' 4 . In 
the rest frame p = it reads: 



k 4 k 



2 i 2 

+ m 



M 2 



M 2 k\ 



<f> E (k4,k) 



a -^L_K E {k^ k; ki, k>)$ E (k' 4 , k') , (12) 



where KE{k 4 , fc; fc 4 , k') — K(ik4,k;ik'^k') and 
$E{ki,k) = $(i/c4,fc). We take sum of ladder, 
eq. 0, and cross ladder, eq. (|10J) . "convert" 
them in Euclidean space and find in this way the 



Table 1 

Coupling constant a = g / '(lGirm 2 ) as a function 
of the binding energy B for ladder kernel with 
H = 0.15 and [i = 0.5 obtained with e = 10~ 6 . 



B 



a((i = 0.15) a(fj, = 0.50) 



0.01 


0.5716 


1.440 


0.10 


1.437 


2.498 


0.20 


2.100 


3.251 


0.50 


3.611 


4.901 


1.00 


5.315 


6.712 



kernel K E = K^' + K y ^ } of the equation |0| 



For the ladder kernel with fi = 0.15 and /j, = 0.5 
and unit constituent mass m = 1 the values of 
binding energy B = 2m — M are displayed in 
Table 2] With all shown digits, they are in full 
agreement with the results of ^21 obtained by 
using the Wick rotation and the method of |13j . 
This demonstrates the validity of our approach. 

The corresponding weight function g for a sys- 
tem with fi = 0.5 and B = 1.0 is plotted in Fig- 
ure |2 Its 7-dependence is not monotonous and 
has a nodal structure; the z-variation is also non 
trivial. We have remarked a strong dependence 
of (7(7,2;) relative to values of the e parameter 
smaller than ~ 10~ 4 , in contrast to high stability 
of corresponding eigenvalues. However, the corre- 
sponding BS amplitude $ and LF wave function 
ip, obtained from 3(7,2) by the integrals (0} and 
0, show the same strong stability as the eigen- 
values. 

The ladder BS amplitude at the rest frame 
p = is shown in Figure |3J On the left, we have 
plotted ^(fco, k) in Minkowski space versus fco- It 
exhibits a singular behaviour due to the poles of 
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the propagators in at fco = ± (efc ± if)) *- e - 
poles moving with fc and M. On the right, we 
have plotted $E(fc4,fc) = <&(fco = «fc4,fc) in Eu- 
clidean space versus k±. Both are calculated by 
eq. J3J with the substitution fco = ?fc4 for the Eu- 
clidean BS amplitude $£(^4, fc). They drastically 
differ from each other though they correspond to 
the same binding energy. The Euclidean solution 
is smooth, in contrast to the Minkowski one. The 
Euclidean amplitude found by direct solution of 
the Wick-rotated equation (|12() is indistinguish- 



able from the one shown at r.h.-side of Figure 
The corresponding LF wave function ip(k±,x) is 
shown in Figure 0| 

The cross ladder effects arc presented in Figure 
where the binding energy Basa function of the 
coupling constant a is shown for exchange masses 
/i = 0.5. Corresponding numerical values - with 
an accuracy of 1% - are given in Table We 
see that for the same kernel - ladder or ladder 
+ cross-ladder - BS and LFD binding energies 
are very close to each other. The BS equation is 
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Figure 4. Wave function 4>(k±,x) for ladder kernel with [i = 0.5 and B = 1.0. On left versus k± for fixed 
values of x and on right versus x for a few fixed values of k± . 

Table 2 

Coupling constant a for given values of the binding energy B and exchanged mass fi = 0.5 calculated 
with BS and LF equations for the ladder (L), ladder +cross-ladder (L+CL) and (in LFD) for the ladder 
+cross-laddcr +stretched-box (L+CL+SB) kernels. 



B 


BS(L) 


BS(L+CL) 


LF(L) 


LFD(L+CL) 


LFD(L+CL+SB) 


0.01 


1.44 


1.21 


1.46 


1.23 


1.21 


0.05 


2.01 


1.62 


2.06 


1.65 


1.62 


0.10 


2.50 


1.93 


2.57 


2.01 


1.97 


0.20 


3.25 


2.42 


3.37 


2.53 


2.47 


0.50 


4.90 


3.47 


5.16 


3.67 


3.61 


1.00 


6.71 


4.56 


7.17 


4.97 


4.91 



slightly more attractive than LFD. At the same 
time, the results for ladder and ladder+cross- 
ladder kernels considerably differ from each other. 
The effect of the cross-ladder is strongly attrac- 
tive. Though the stretched box graphs arc in- 
cluded, their contribution to the binding energy is 
smaller than 2% and also attractive. This agrees 
with the direct calculation of the stretched box 
contribution to the kernel performed in |14j . The 
solution of the Wick rotated equation 112|) for 
L+CL kernel was found in |15| . Within a numer- 
ical accuracy of 1%, it coincides with the corre- 
sponding solution given by our method, i.e., with 
the column BS(L+CL) of the Tabled This coin- 
cidence proves the possibility to solve BS equation 
with cross ladder kernel by using the Wick rota- 



tion. It also confirms the validity of our method. 
We emphasize again the fact that the Wick ro- 
tated equation allows to find the binding energy, 
the same as the one determined by the Minkowski 
BS equation, but not the Minkowski BS ampli- 
tude. 

6. CONCLUSION 

We have developed a new method for solv- 
ing the Bcthc-Salpctcr equation in the Minkowski 
space, i.e. without making use of the Wick rota- 
tion. This method is based on an integral trans- 
form of the original BS equation which removes 
its singularities. 

Our equation is more easily formulated in terms 



r 




a 

Figure 5. Binding energy B versus coupling 
constant a for BS and LFD equations with 
ladder (L), laddcr+cross-laddcr (L+CL), and 
ladder+cross-ladder+stretched box (L+CL+SB) 
kernels. The exchange mass is fi = 0.5. 

of the weight function of the Nakanishi integral 
representation [|J. This allows a straightforward 
reconstruction of the original Bcthe-Salpeter am- 
plitude - both in Minkowski and Euclidean space 
- as well as the Light-Front wave function. 

For zero-mass ladder exchange, we reproduce 
analytically the Wick-Cutkosky model. For mas- 
sive ladder exchange, the binding energies are 
found numerically and are in full agreement with 
preceding results obtained in the Euclidean space. 

For a given kernel, Bethe-Salpeter and Light- 
Front approaches give very close results, though 
the first one is always slightly more attractive. 
The cross-ladder contribution is strongly attrac- 
tive in both models. The strctched-box contri- 
bution - evaluated separately only in the Light- 
Front framework - is attractive but small. 

Binding energies obtained by solving Bethe- 
Salpeter equation in Minkowski space with ladder 
plus cross-ladder kernel are also in full agreement 
with the ones obtained in Euclidean space |15j . 

We are aware about only one work |16| where 
the separated effect of the cross-ladder diagrams 
in the Bethe-Salpeter equation has been esti- 
mated. The authors used an approximate disper- 
sion relation to obtain the corresponding kernel. 
Our results are found to be three times smaller 



than those given in this reference. 

Our method for solving the Bethe-Salpeter 
equation in the Minkowski space can be gener- 
alized to non-zero angular momentum and, pre- 
sumably, to the fermion case. 
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